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Abstract 

Consider a non-symmetric generalized difTusion X{-) in determined by tiie difTeren- 
tial operator A{x) = — . diaij{x)dj +'^- hi{x)di. In tliis paper the diffusion process 
is approximated by Markov jump processes Xn(-), in homogeneous and isotropic grids 
G„ C R , which converge in distribution to the diffusion X[-). The generators of Xn(-) 
are constructed exphcitly. Due to the homogeneity and isotropy of grids, the proposed 
method for d > 3 can be apphed to processes for which the diffusion tensor {aij{x)}ff 
fulfills an additional condition. The proposed construction offers a simple method for 
simulation of sample paths of non-symmetric generalized diffusion. Simulations are 
carried out in terms of jump processes Xn[-). For d = 2 the construction can be easily 
implemented into a computer code. 

AMS subject classification: 60H35, 60J22, 60J27, 60J60, 65C30 
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1 INTRODUCTION 

A symmetric tensor valued function x a{x) = {aij{x)}f'l which is measurable, bounded 
and strictly positive definite on R"^ defines a second order differential operator in divergence 
form on K'', Ao{x) = — J^ij diaij{x)dj. Each Ao(^c) determines a symmetric diffusion X{-) 
in . In 'SZ the process X{-) is approximated by a sequence of Markov jump processes 
Xn,{-) (MJP) in grids, i.e. processes in continuous time and discrete state spaces. Coordinate 
systems in in |SZ| depend locally on the structure of the tensor valued function a. The object 
of the present analysis is a non-symmetric diffusion determined by A{x) — Aq{x) -|- B{x), 
where B{x) = ^j6,i(x)i9i, and a construction of its approximations by MJPs in lattices. 
For d — 2 the method is valid for any a, and for d > 2 the method is valid for tensor valued 
functions a satisfying an additional constraint. In this way we offer an efficient method 
for simulation of sample paths of a non-symmetric generalized diffusion using MJPs. Each 
Xn{-) can be simulated by well-known methods. 

In the case of classical diffusion determined by an elliptic operator of the form A{x) = 
— ^j^- aijix)didj + bi{x)di, where , bi are Holder continuous on R"^, approximations by 
MJPs can be efficiently used to simulate the first exit from a bounded set of R"*. However, 
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such an approach is one of several existing possibihties, and the motivation for the construc- 
tion in terms of MJPs happens to be of lesser importance, since the process has a repre- 
sentation in terms of SDE which can be simulated straightforwardly. For one-dimensional 
generalized diffusion, there exist representations in terms of SDE as described by Etore |Et| 
and Lejay & Martinez |LM] . so that such representations can be used for simulation. In the 
case of a process defined by a differential operator in divergence form on M'', d > 2, there 
is no such natural representation, so approximations by MJPs are essentially the only tool 
available for simulations. 

A class of convergent approximations Xn{-) in [SZ] is constructed by using discretizations 
of the corresponding Dirichlct form {v,u) i— >■ a(v,u). The functions diV,diU in a{v,u) are 
approximated by forward difference operators in local basis which generally varies. This 
approach is anticipated in ^MW) without any remarks on the convergence of the constructed 
MJPs Xn{-)- In our approach, the MJPs Xn{-) are constructed in terms of generators 
A„(gen) on the grids, 

d 

Gn = {hn J2 ^'^^ • ^' ^ 2-", (1) 

with a fixed basis {e^}^. The index set of grid-knots is denoted by /„. In order to sim- 
plify expressions, we often write h instead of The generators A„(gen) are constructed 
explicitly from a general principle which is not directly related to forward difference oper- 
ators. Then the discretizations h'^an{v,u) of the original Dirichlet form are associated to 
the constructed generators yl„(gen). It turns out that h'^an{v, u) cannot be simply obtained 
from discretizations of original form by using forward/backward difference operators. The 
obtained class of A„(gen) is not included among generators constructed in |SZ| . 

Advantages and drawbacks of the present construction can be briefly described as follows. 

Advantages: The convergence of MJPs is proved for a non-symmetric generalized diffu- 
sion. The generators A„(gen) are explicitly given in terms of values of functions aij, hi. For 
the case of d = 2, the matrix entries of A„(gen) can be easily implemented into a computer 
code because rotations of coordinates are avoided. For d = 2, the construction holds for a 
general matrix valued function a(-) on R^. 

Disadvantages: For d > 3, the proposed construction is not valid for all a(-) on M''. The 
restriction is defined by an auxiliary matrix valued function a(-), with the following matrix 
entries: 

du = an, dij = -\aij\, i ^ j. (2) 

The here proposed construction of A„(gen) is valid only if a(-) is strictly positive definite on 
W^. For d > 3, there are simple examples of pairs a, a, where a is definite and d is indefinite. 
We can say that ([2]) is valid if the off-diagonal entries of a(-) are small in comparison with 
the diagonal ones. When d{x) becomes indefinite, it is necessary to apply a local rotation of 
coordinates, that ensures a diminishing of the off-diagonal entries, thus ensuring the positive 
definiteness of d{x) in the new coordinate system. 

Now we can describe basic steps in the construction and the proofs. Let U{-) be the 
strongly continuous semigroup in the Banach space C'(R'^) (continuous functions vanishing 
at infinity) which is associated with the diffusion process X{-) |EK] . For each n the space 
of discretizations of (^(R^) in terms of grid-functions on Gn is denoted by /oo(G„). The 
generators A„(gen) determine semigroups [/„(•) in ?oo(G„). There exist continuous mappings 
*n : /oo(G„) ^ CiW^) and : CiW^) k> /oo(G'n) with the following properties: |1$„|| = 
1, < I and $,7i$„ = / on i^{G„). For / e C(R''), the grid-function f„ === <i>-\f is 

called the discretization of /. Now, for each t > 0, the grid-functions u„(t) — ^^^U{t)f and 
Un{t)in can be compared. We need the following relation: 

lim sup{|C/„(t)f„-$-i;7(t)/|oo : te [0,1]} - 0, (3) 

where | • |oo is the Zoo(G„)-norm. By Theorem 2. II of |EK| . this relation is a sufficient 
condition for the convergence in distribution of MJPs Xn{-) to the diffusion process X{-). 
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Definitions of various objects are given in Section [2] In Section [3l a class of generators 
^„(gen) is constructed explicitly. Grid- functions u„(t) = exp(v4„(gen)t)f„ and their images 
u{n,t) — $„u„(t) in C'(R'') are studied in Section SI The convergence u{n,t) U{t)f 
in C'(M'^) is proved initially for the case of smooth coefhcients Uij^bi, and then for general 
coefficients a^, 6^. A numerical example is given in the last section. 

2 PRELIMINARIES 

The Euclidean norm in is denoted by | • |. All the open subsets of W^, considered in this 
work, are bounded and connected open sets with Lipshitz boundary |Mal ISe] . We call a 
subset of this kind a Lipshitz domain, denote it by D, and its boundary by dD. 

The Banach spaces of functions C'^''^R'^),C^/\r'^) = Co(M'') n C^'''>{R'^) are defined as 
usual, Co(ffi'^) being the linear space of continuous functions with compact support. Their 
norms are denoted by || • The closure of functions in C(M.'^) with compact support 

determines the Banach space C{M.'^). The Holder space of parameter k+a, k G No, a G (0, 1), 
is denoted by C^'^+^^R'^) and defined as the completion of C^^^^W^) in the norm: 

Nl'-"^ = NI^) + sup{ '"^-(-+g-"^"(-)^ .eM^o<N<i}. 

The Lp-spaces as well as Sobolev -spaces are defined in a standard way }Ma[ ISe] . Their 
norms are denoted by || • \\p and || • ||p.i, respectively. For each p, 1 < p < oo, the norm of 

W^iR') is defined by \\u\\,,^ = {Ml + \\\7u\\lf\ where || Vu||p = (E.ti \\d,u\\lf'. 
We say that a function / on R'' belongs to a class C(''+") on R'^ if < oo. 

In this article we consider a 2"'^-order elliptic operator on R'^, 

d d 

A{x) = - J2 9^a^J {x)d, + J2 h (^) ^3 > (4) 
for which the coefficients fulfill the following: 

Assumption 2.1 The functions aij = aji, bi,i,j — l,2,...,(i, are measurable on R"^ and 
have the following properties: 

a) There exist positive numbers M, M, < M < M, such that the strict ellipticity is valid: 

d 

M\z\^< a,j{x)z,Zj <lA\z\^, x,zgR'^. (5) 

ij = l 

b) The functions bi are bounded on R"^ . 

In our analysis we regularly use the notation Ao{x) = — J2ij=i ' ^i^) = 
Y^j^ibj{x)dj and A{x) = Ao{x) + B{x). The following real bilinear form on W^{D) x 
W}{D), 1/g = 1, defined by: 




is associated with the differential operator A{x). 

A basic result towards a proof of ([3]) is the following theorem [St] : 

THEOREM 2.1 Let the differential operator A{x) be defined by ^ and Assumvtion \2. 1\ 
Then —A{x) has the closure in CiR!^) and generates a Feller semigroup in CiR!^). 
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The semigroup of this theorem has the restriction to a closed subspace of C''-"^(M'') = 
C'(M'^) n C''°'\W^) which is also a strongly continuous semigroup. 

COROLLARY 2.1 Let A{x) be defined as in Theorem \2.1\ There exists an a & (0, 1) and 
(T > that depend only on M_,M and ||b||oo, such that the following two assertions are valid: 

(i) The operators U{t) in C^"^(W^) are bounded uniformly on segments o/[0, oo). 

(ii) There exists a closed subspace i^'^")(R'^) C C''-"-' (M'') such that the closure of —A{x) in 
p{o) (R"*) generates a strongly continuous semigroup U{-) in F^'^^iW^), with \\U{t)\\\S,^ < 
exp(crt). 

Proof: The Feller semigroup of Theorem 12.11 can be represented as t exp(t)F(t), 
where V{-) is the strongly continuous semigroup generated by the closure of — (/ + A{x)) 
in C'(IR.''). Due to Theorem 12.11 the differential operator / + A{x) has the bounded inverse 
(/ + in C{M!^). For any / e C{R'^) the continuous function {t,x) ^ u{t,x), defined by 
u{t) = V{t)f, is a solution to the following IVP for PDF 

dtu{t,x) + {I + A{x))u{t,x) = 0, (7) 

with the initial condition 7i(0, x) ~ f{x). In addition, ||u(t)||oo < WfWoo for all t>Q. It turns 
out that a solution to (O, for each t > 0, is an element of C*^"-' (R'^). For an initial condition 
/ S C^"-* (R"*) we have a stronger result, a solution to ([7]) is an element of C^"-' (R^) uniformly 
with respect to t G [0,1]. This follows from the following arguments. Let us consider the 
balls Bi{x) C M'', centered at a; e M'' and having the radius equal to 1, and let us consider 
the restrictions u(t)\Bi{x) for any a; g M'* and any t G [0, 1]. According to the basic theorem 
of Section 10, Chapter 3 of |LSU| . there exists a € (0, 1) depending on M_,AI, ||b||oo and /3 
depending on M, M, |ib||oo and ||/||^^ such that a solution m to ([7]) on Bi{x) is an element 
of C(")(Bi(a;)), and 

sup ||^.(^)|Sl(rr)|jM + sup sup ~ "f"' < /3. (8) 

te[o,i] t,se[o,i] yeBi{x) 

This inequality implies that the operators V{t) ,t € [0,1], are linear in C*^"^ (R'') , with a norm 
which is uniformly bounded with respect to t S [0, 1]. Hence, the operators V{t) define a 
semigroup of bounded operators in C'^°'\R'^). This proves the assertion (i) of the corollary. 

For any g e C'(R'^) the function / = {I + A)^'^g belongs to the domain D{A) of the 
closure of A{x) in C(R'^), and V{t)f = (I + A)-^V{t)g. The linear space spanned by 
f = (I + A)~^g, g e C'(R''), is denoted by E(R'^). Let us assume now that we can prove 
f = (I + A)-\j' e C'(")(R'') for any g e CiW^), i.e. < P\\g\\oo, where /3 is a 

number depending on M, M and ||6||oo- This would imply that i?(R'') consists of elements 
in C'^"^(]R'') and its completion in the || • ||(")-norm is a closed space 
In addition, the above assumption would imply the following inequality: 

\\iVit)-I)f\\^-^ < (3UVit)~I)g\U (9) 

i.e. the continuity of the function 1 1-> V{t) on a dense linear space bmce 
the operators V{t) are bounded in C^°'^{M.'^) uniformly with respect to t G [0, 1], the obtained 
inequality would imply that the operators V{-) define a strongly continuous semigroup in 
the Banach space F(")(R''). Then the semigroup U{-) of (ii) is also strongly continuous on 
i^(")(R'') and the existence of cr > in (ii) follows from the theory of strongly continuous 
semigroups. 

Therefore, in order to prove (ii) it remains to show that (/ + A)^^ maps C'(R'*) into 
C'(")(M''). For this we consider the elliptic problem: 

{I + A{x))u{x) = f{x), xeR"^, (10) 
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where / £ C'(")(R'^). The existence of sohition u = {I + A)-^ f in C{R'^) is granted by 
Theorem 12. 11 and in addition we have ||u||oo < ||/||oo- In order to prove that u G C'(")(R''), 
it is sufficient to verify that u\Bi{x) S C'^"-'(i3i(a;)) for any x and a fixed a. To obtain this, 
we apply the theorem in Section 14, Chapter 3 of |LUj . The fohowing estimate is vahd: 



< /3 11/11 



oo ■ 



where both, a g (0, 1) and (3, depend on M, M, \\b\\ 

The parameters a S (0, 1) for the parabohc (0 and elliptic problems (ITUl) are not neces- 
sarily equal. We choose the minimum to obtain ([9]), and consequently the assertion (ii). QED 

The assertion (i) of Corollary 12.11 is used in the proof of Theorem 14. 1[ which establishes 
the basic inequality ([3]). 

The linear space of grid-functions on G„ is denoted by l{Gn)- Elements of /(G„) are 
also called columns. Columns are denoted by u, v etc, while their entries are denoted by 
UkjVi etc. Thus a column u„ has entries (u„)fc,fe G /„. Columns with a finite number of 
components span a linear space lo{Gn)- The completion of lo{Gn) m the | • |oo-norm is 
denoted by l{Gn). The corresponding Z^-spaces are denoted by lp{Gn) and their norms by 
I • \p. For p < oo this norm is defined by |u|p = [X^fce/ ^-^d for p = oo by 

|u|oo = sup{|ufc| : k e Z''}. The scalar product in hiGn) is denoted by ( -j- ) and sometimes 

by(-l-)- 

Let / e C(M'*) and the column f„ £ l{Gn) be defined by its components (f„)fc ~ f{hk), 
where h = hn — 2^". Then the mapping C{W^) 3 / H> f „ G l{Gn) is called the discretization 
of /. 

The shift operator Z{x),x G M*^', acting on functions / : M*^' M- M, is defined by 
{^Z{z)f^{x) = f(x + z). Similarly we define the discretized shift operator. The shift 
operator Zn{r,i) by r units in the direction is defined by (Z„(r, j)u„)^ = (u„)i, where 

I = k + rSi. The partial derivatives of u G C'-^^(M'') with respect to a grid step h are 
discretized by forward/backward finite difference operators in the usual way, 

a^{th)u{x) = ^{u{x + thei) - u{x)), ^^^d^^^Q^ 
Bi{th)u{x) — ji^iu{x) — u{x — thei)^, 

Let r G Z \ {0}. Discretizations of the functions diU on G„, denoted by C/i(r)u„, Vi(r)u„, 
are defined by: 

{Ui{r)Un)^ = Bt{rh) u{Xm.), {Vi{r)Un)^ = Bi{rh) u{xm), 
where a;™ G G„. Then 



Udr) - {rh)-^Zn{r,z) - I), 

V,ir) = (r/i)-i(/-Z„(-r,z)) = U,i-r) = -U,{rf 



(12) 



Therefore we have Ui{—r) = Ui{r) Z„(— r, i) = Zn{—r,i) Ui{r), and similarly for Vi{-). We 
use the abbreviations Ui — Ui{l),Vi = Vi{l). 

A matrix A„ on G„ is called a matrix of positive type if the diagonal entries of A„ are 
positive, off-diagonal entries are non-positive and the row sums are non- negative. If An{gen) 
is the generator of a MJP in G„, then —An{gen) is a matrix of positive type. 

3 CONSTRUCTION OF GENERATORS OF MJPs 

To discretize A{x) means to associate to A{x) a sequence of matrices A„ on G„, n £ N, with 
the following properties: 

a{v,u) = lim /i'^(v„|A„u„), v, u G G^'^ (R"^). 
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The terminology "discretizations" of A{x) instead of approximations of A{x) appears to 
be more suitable at the beginning of the construction, since the convergence analysis is 
postponed to Section IH 

We wish to emphasize that discretizations An of the differential operator Ao{x) are 
derived from a general principle, similar to the one exploited in |LR2| . This method is not 
based on finite difference formulas. Nevertheless, bilinear forms need to be associated to An 
so that An can be derived from the corresponding variational equalities. The constructed 
bilinear forms can be considered as the discretizations of the original form ([6]) . These forms 
are basic objects in our proof of convergence of discretizations. 

Discretizations to be considered in this section are possible if certain conditions on 
are fulfilled. The required conditions are stronger than those given in Assumption 12.11 By 
relaxing them gradually as n — >■ oo we obtain discretizations for a general A{x) given in 
Assumption 12.11 

To each pair t; e Gn,P& N*^, there is associated a rectangle C„(p, v) =Ylf^i[vi,Vi + hpi) 
with the " lower left" vertex v and the edge of size hpi in the i-th coordinate direction. These 
rectangles define a partition of M''. Apart from these rectangles, we will need the closed 
rectangles, 

d 

which are centered at the grid-knots v. Evidently, Sn{p, v) is the union of closures of those 
rectangles Cn{p,x) which contain the grid-knot v. 

A discretization An is defined in terms of its matrix entries {An)ki, where hk, hi G G„. 
For a fixed x = hk G Gn the set of all the grid-knots y = hi such that {An)ki ^ is denoted 
by Af{x) and called the numerical neighborhood of An at x G Gn- The set Af{x) contains 
always a subset consisting of x and 2d elements x ± hei, i = 1,2, ... ,d. Additional elements 
ofAf{x) are possible as we shall see, depending on the sign of ,i j. In terms of the MJP 
Xn{-), the set Af{x) consists of the states of possible jumps from the state x. Let us point 
out that the sets J\f{x) vary with n, that is for two grids G„, Gm, n ^ m and x e Gn n G™, 
the corresponding numerical neighborhoods Af{x) are different. 

General setup 

In order to give a comprehensive insight into the proposed construction of the discretizations, 
it is convenient to initially consider a differential operator Aq — o-ijdidj having a 

constant diffusion tensor a = {ajjlii- In this case, the matrices An need to have a property 
which is called the consistency. Let x M- p{x) — po + J^iPi^i + Pij^i^j^ Po,Pi,Pij G 
be a second degree polynomial in arguments Xi , and let p„ be its discretization on grid G„ . 
Then the consistency holds if the following identities are valid: 

d{A,n,x)p{x) := A{x)p{x) - {AnPn) ^ = 0, x = hk G Gn- (14) 

These consistency conditions are sufficient for proving the convergence of L/„(t)f„ to $„J7(t)/ 
in the Banach space of continuous functions on as specified in Actually, the first 
step of construction of matrices An begins with a search for those matrices An which are 
simultaneously of positive type and fulfill the conditions S{A,n,x)p{x) — on G„ [LRl . 
We request that the discretizations An have the following properties: 

a) The numerical neighborhoods A/'(a;) C G„ resemble each other, that isAf{x) = x+Af{0). 

b) The numerical neighborhoods J\f{0) C G„, n G N, resemble each other, i.e. ii y = hnl G 

A/'(0) C Gn then y' = /i,„Z G U{0) C G^., n,m G N,n m. 

c) The matrices An are symmetric. 

d) The matrices An are consistent discretizations of Aq. 

Now we have the following result. 
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LEMMA 3.1 If a is positive definite then there exist matrices An of positive type fulfilling 
the conditions a)-d). The non-trivial entries of An are defined in terms of d natural numbers 
ri , r2 , . . . , G N hy the following expressions: 

(^n)fcfc-j_e^ = -773- - Hm^i 7^\o-im\ , «=l,2,...,d, (15) 

i^^) kk±z{i,j,±) ^ ^ TT^TITJ l^-yl' j = 1, 2, . . . , d, 

where in the last line, z{i,j, +) — riBi+rjEj corresponds to the case Oij > 0, and z{i,j, —) — 
TiEi — rjBj corresponds to the case Oij < 0. 

Proof: For a matrix An defined by (|15p the conditions a)- c) are obviously satisfied. 
(Recall that Af{hk) consists of those grid-knots hi E JV{hk) for which (A„)^j are non- 
trivial). It remains to prove the condition d), and that An is of positive type if a is a 
positive definite matrix. 

The condition d) is equivalent to the following property. Let x i-^ p{x) — Pij^i^i + 
'^iPiXi -|-po be a second degree polynomial. Then d) is valid if and only if 

hke^f{o) i<j 

For An defined by (|15p this identity can be easily verified by using the monomials p{x) = x^Xj 
with i ^ j and p{x) = a;|, for aU i,j — 1,2, ... ,d. 

In the last step of the proof we show that there exists a sequence ri,r2, ■ . . ,rd £ N 
for which the brackets in the second line of (|15p are positive, thus ensuring the matrix An 
to be of positive type. For this purpose we assume that the matrix a is positive definite. 
Let us consider the eigenvalue problem aw = Xw, and let us also assume that the ma- 
trix a is irreducible. For /i > sufficiently large the irreducible matrix fil + d has the 
inverse (/i/ + a)^^ with positive entries, so the Perron-Frobenius theorem can be applied to 
(/z/ + a)~^. Thus the eigenvector corresponding to its maximal eigenvalue is positive. This 
result can be also formulated in terms of the problem aw = Xiw for the minimal eigenvalue 
Ai of a. We have Aiii; > and consequently dw > 0. This inequality can be rewritten as 
O'ii ~ 'l2m^ii'li / 1m)\o-im\ > whcrc Qi = . One can find rational approximations ri/rm 
of Qi/qm which also fulfil the obtained inequalities. If d is not irreducible, it can be rewritten 
in a block diagonal matrix form with irreducible blocks. The previous construction can be 
applied to each block. QED 



In the case of d = 2, either both a,d are positive definite or neither is. For d — 3, 
there are symmetric positive definite matrices a for which d are indefinite. For instance, 
the symmetric matrix a of order d — 3 defined by an = 1, ai2 = 023 = — 1/\/2 has positive 
eigenvalues for the case of ai3 > and a negative eigenvalue for the case of ai3 < 0. 

For d = 2 and ri = 3,r2 — 1, two possible numerical neighborhoods Af{x) are illustrated 
in Figure [T] 




ai2 > ai2 < 



7-1=3, r2 = 1 
Figure 1 

Numerical neighborhoods 
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If An is a matrix from Lemma 13. 1[ then An{gen) = is a generator of a MJP in 

Gn- Let us consider now functions on R'' defining a diffusion tensor at each a; € R''. 
Let us assume in addition that the matrix a{x) — {aij{x)}ff is positive definite at each 
X S ISJ^. One could replace the numbers a^- of ([T5|) by the numbers aij{hk). The result- 
ing An{gen) — —An would again be a generator of a MJP in Gn- However, the matrices 
An thus obtained are discretizations of the differential operator A{x) = aij{x)didj, 
as will be seen in Section SI and not of A{x) = Here we present a 

method of construction of discretizations of a differential operator in divergence form, 
A(x) = — ^j^j diaij{x)dj, resulting in matrices of positive type, with a structure similar 
to (jlSp . The operator in divergence form, A{x), naturally corresponds to the bilinear form 
a{v,u) = J aij{x)div{x)dju{x)dx. Consequently, its discretizations An will correspond 
to a sequence of discretized forms a„(w, u) on grids G„. 

Construction for d = 2 

Initially one considers the constant coefficients and constructs the forms an{v,u) — 
(v„|A„u„) for which the matrices An coincide with (IT5|) . Then the obtained expression 
of an{v,u) is generalized to the case of non-constant coefficients aij. The corresponding 
matrices An are obtained from the variational method in the standard way. 

For the case of constant coefficients a^j, the bilinear form /i^a„(u, u) that discretizes 
the form a{v,u) = ciijdivdjU, will be a second order polynomial in the quantities 

Bi{rih)v{x), Bj{rjh)u{x) with a certain choice of ri,rj € N and x G G„. In order to write 
down the form as simply as possible, we use the following abbreviations: 

Ukihfi) = 3i{rih)u{hk) = {rih}^^ [u{hk + Vihei) — u{hk)], 
Ukii,r.i) = 3i{rih)u{hk) = {rihy- [u{hk) - u{hk - rihe^)]. 

For ai2 < the form is defined by: 

( ' 

al Hv,u) = ^ ^ aiiVk{i,l)uk(i,l) + ^ a^jVk{i,ri)uk{j,rj) 

k \i=l i^j 

+ a-ij— [vk{iA)uk{iA) ~ Vk{i,ri)uk{i,ri)\ 

If the last summand of ^TE\\ were omitted, then the resulting An would have neighborhoods 
N{x) containing seven grid-knots x,x ± he^, x ± h{riei — r2e2) as illustrated in Figure [TJ 
as well as the following four additional grid-knots, x ± hviBi. The second line of ([TC]) causes 
cancellation of those matrix entries which would correspond to four additional grid-knots. 

For ai2 > the form is obtained from the previous one by changing Vk{i, •), Uk{i, ■), i — 2, 
into t}fc (?,•), •), i = 2, respectively. 

For the case of functions x i— >• aij{x) one naturally starts from the just obtained expres- 
sion for ait\v^ u), changing the numbers aij into the numbers aij{hk). In fact, the numbers 
aij{hk + t), where t are appropriately selected element of R^, are acceptable. The fastest 
convergence of h\an{v,u) — >■ a{v,u) is a criterion which helps us to choose the vectors t. It 
turns out that the best choice is t(r) = {h/2){riei + ^262) |TS| . i.e. the values for which 
hk + t{r) are the mid-point of the rectangle Cn{r, hk), with the lower left vertex at hk and 
the upper right vertex at hk + riBi -1-^262. Let us remark that the finite difference operators 
3i{rih)u{hk) are defined in terms of function values at the vertices of Cn{r,hk). Thus, if 
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cii2{x) < 0,3? G M^, we get 



T,t^j aij{x + t{r))vk{i,n) Uk{j, rj) 
+ ^ aij{x + t{r)) — [vk{i, l)uk{i, 1) - Vk{i,r,) Uk{i,ri)]) , 

where r = (ri,r2),l = (Ijl)- By using the variational method, one obtains the entries 
(^n)fei,fe,' € in- At the present step of construction it is not necessary to write down all 
the entries. In order to describe the influence of the parameters r\,ri on the structure of 
entries we consider one group of entries: 



1 



The matrices An are of positive type iff the bracket has positive sign for each x € Gn- This 
is a condition on the functions aij , which is imphed by a particular choice of ri , r2 G N. 

For ai2 > the form is obtained from the constructed one by changing the following quan- 
tities. The finite differences Wfe(i, •), Mfc(«, •), i = 2, should be changed into t;fc(i, •), Ufe(i, = 
2, as in the case of constant coefficients, and t(rn) should be changed into s(m) — {h/2){miei~ 
771262). With these changes, the obtained An are symmetric matrices, possibly of positive 
type. We say "possibly of positive type" since this depends on the choice of ri , r2 . 

In the general case, the sign of ai2 is not constant on R^. Therefore, we partition 
into two subsets {x £ : ai2{x) < 0} and {x £ : 012(0;) > 0}. Then each of these 
sets has to be partitioned further, where each of the partitioned classes is characterized by 
a pair ri,r2, so that the resulting entries {An)kk±ei have negative values. The construction 
is carried out for a class of functions Uij with moderate discontinuities. 

Assumption 3.1 Let there exist a finite index set L, a partition R^ = U/g^D/ and a 
diffusion tensor a = {aij}i\ on R^ satisfying the strict ellipticity conditions ^ and the 
following additional discretization conditions: 

a) There exist q G N^, no G N and the corresponding ha = 2~^° so that each set Di is a 

connected union of cubes of the form Y\^^i{xi,Xi +qiho). The matrix-valued function 
X n- a{x) is continuous on cl{Di) and either 012 > or ai2 < on Di . 

b) To each Di there is associated a parameter r{l) G N^, such that the following inequality 

is valid: 

uj(a) — inf min min inf < inf duiz) 

n l£C i=l,2 L zeS„(r(i),x)nA 

nil) \ (1"^) 

^-ryr SUp |a„„(,;) : X <£ Gn\ > 0. 

where m(i) — 3 ~ i and the rectangles Snii'{l),x) are defined by 113\) . 

Condition b) is crucial in our construction of discretizations An with the structure of 
matrices of positive type. 

The set L of Assumption 13 . 1 1 is partitioned into the subsets L^, where I G C_ means that 
ai2 < on Di, and I E means that ai2 > on Dj. It is convenient to use a representation 
an{v, u) = ail \v,u) + ai^\v, u), where the form ol \v, u) contains the sums over the grid- 
knots in Di, I G £j- and a^rt\v, u) contains the sums over the grid-knots in Di, I G Let 
us define 

f(±+)(^) ^ |(± ri(Oei + r2(0e2) G5„(r,0), 

t(±-)(r) = |(± ri(Oei - r2(0e2) G 5„(r,0). ^ ^' 
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where r — (ri,r2). The form ai ■*(•,•) is defined by: 



ai \v,u) 



au{l,x + t(++\l)){d,{h)v){x){a,{h)u){x) 

i=l 

a„- il,x + t(++) (r)) ( a (a;) ( 3 ^ (r, (a 



... { E 
+ E 

+ E 

- {B,{r,{l)h)v){x){B,{r,{l)h)u){x)\ )} 



a 



(/, a; + (r)) 



a,{h)v){x){a,{h)u){ 



X 



(19) 

where, as aheady noted, r = (ri,r2), so that 1 = (1,1). The form ai^\v,u) is obtained 
from ail \v,u) formaUy by replacing with £+, then 3{mh)f with 3{rah)f for each 
'T^ G {li^'z} and / = and with t^^ Observe that the forms a^r^\v^u) are the 

second degree polynomials in the quantities 3i{qh), 3i{qh) with q e {l,ri,r2}. 

We say that x £ Gn is an internal grid-knot if M{x) C G„ n Di for some All the 
other grid-knots are called boundary grid-knots. For an internal grid-knot x the acceptable 
expressions for {An)ki follow directly from the definition of corresponding discrete forms 
V, u i-^- (v|A„u). For a boundary grid-knot x = hk the obtained expressions are complex, 
and the calculated (j4„)fcj could break down the structure of matrices of positive type. 
Therefore, one seeks simpler procedures for constructing entries (A„)fej for boundary grid- 
knots X — hk. The results of such a construction must be matrices An which determine 
MJPs and the convergence ^ of MJPs determined by A„ should also be ensured. 

In order to write down the entries of A„ corresponding to internal grid- knots, we use 
and the following abbreviations: 



(0 



w 

(a/3)/ \ 



,(++) 



(r) 
(r) 
(r) 



ri(/)ei ± r2(/)e2 G /„, 
a,,(a; + t("«(r)), a,/3e {+,-}, 
ai2ix + t^++\r)-hei), 
ai2ix + t^++Hr) -he2), 



^(++)(r), a(r-)(r) = a^->ir). 



When we apply variational method to the form ^ defined by (flQl) and the corresponding 
form ai^\ we obtain entries of A„. Thus, for an internal grid-knot hk the nontrivial off- 
diagonal entries of A„ are: 



{Ar, 



{Ar, 



' fcfcibei 



1 
1 



'■11 



(1) 



^22 



(±-) 



(1) 



(1) 



r2(l 

_ nil 

r2{l) 
r2il) |, 

nil) 



a 



(±+) 

12 



(r-)l 



|«i2 

j( + ±) 
'■12 



\r)\ 



for ai2 < 0, on D;, 
for ai2 > 0, on D;, 



(20) 



(r)| 



The entries corresponding to the grid-knots in the plane spanned by ei , 62 have the structure: 

1 



(^")fcfc-«,(-)(/) 

("^")fcfe + TO( + )(/) 

("^")fcfe-tc(+)(i) 



/iVi(/)r2(/) 
1 

/iVi(j)r2(0 
/i2ri(/)r2(/) 
/»2ri(0r2(0 



*12 



-^12 
-,(++) 



'{r)\ 
{r)\ 
{r)\ 
{r)\ 



for ai2 < on Z?; 



for ai2 > on £)/ 



(21) 
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If the quantities dl°'^\r),a, l3 e {+, — } in ([20)) are replaced with a-J'^^(r), the conver- 
gence is stih preserved. However, the quantities a\°'^\r),d\°'^\r) should not be replaced 
with aij{hk) since the resulting {An)ki would be discretizations of — cnjdidj and not of 
— diGijdj. This assertion can be easily proved for the case of dimension d — 1 and the 
diffusion tensor a, that is twice continuously differentiable. We intend to compare the ex- 
pressions A{x)u{x) = —{a{x)u{xyy and A'^^{x)u{x) = —a{x)u"{x) and their discretizations 
on grids G„ = {hk : fc e Z} C M. The discretizations of A'^" {x)u{x) are given by Lemma 



- -h-^a{hk), (A-)^^ = 2h-^a{hk), 

and consequently: 

2u{hk) - u{hk + h) - u{hk - h) 



((A^^)„u„)^ = aihk)- 



For the discretizations of A{x)u{x) we consider (pHl) . Let us represent a{x ± h/2) by its 
Taylor polynomial of the second degree, a{x ± h/2) — a{x) ± {h/2)a'(x) + {h/2)'^a" (x) + 
r(±/i, x), where the remainder r(±/i, x) has the property lim^_^o h^'^r{±h, x) — 0. Therefore 
we get: 

where lim;j_>o 2;) = 0. Hence: 

. / h'^a"(hk)\ , ^ 2u(hk) - u(hk + h) - u(hk - h) 

+ ,^^f^k f^hk + h)-^uihk-h) ^ 

From the obtained expressions we have: 

(A„u„), - + = (Snu„)^, +5(/i,a;), 

where (i3„u„)j^ are the discretizations of ~a'{x)u'{x). In other words, we have obtained 
discretizations of the expression A{x) — ^'^^(2:) — a' {x){d/dx) as we should have. 

Now we can describe the construction which gives a satisfactory result for both types of 
grid-knots, internal and boundary ones. For each a; = /ifc e cls(Z)/) the entries {An)ki are 
constructed by the rules (PO)) . (PT|) . If M{x) C cls(Z?;) there is nothing more to adjust. If 
N{x) n cls(£'/) ^ N[x) then the quantities a"^^\r), d^^^\r), where a,l3 <E {+, — }, should 
be replaced by zeros in all the expressions (|20]) . ([2T|). Let us point out that otherwise the 
entries (|20l) for the case of N{x) n cls(D/) ^ A/'(a;) could have a wrong sign, i.e. it could 
happen that {An)ki > for some k ^ I. In order to avoid such undesired features, one 
therefore omits the terms of entries in the expressions (|20l) , (pTjl , which are proportional to 
|ai2|. This adjustment procedure is equivalent to the assumption that the function 012 is 
zero in a neighborhood of the boundary F = UidDi. 

This determines the rules of construction of discretizations of a generalized diffusion for 
which the diffusion tensor satisfies Assumption 13.11 The above described construction of 
entries {An)ki at boundary grid-knots can be justified as follows. 

A numerical neighborhood Af{x), where x G Di Gn,l G £ is contained in the closed 
rectangle 

d 

Sn{r{l),x) = ]J [.Xi - /i„ri(Z) 

.(0]- (22) 

The union of all such rectangles centered at boundary grid-knots x gT n Gn is denoted by 
S'„(F). This is a closed set covering F. Because of Gn C Gn+i we have S'„+i(F) C 5'„(F) 
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and the identity F = n„S'„(r). For each n e N, we approximate the original diffusion tensor 
a by the tensor a*^"-' defined as follows: 

a,j{x) for xeM.^\Sn{r), ije{l,2}, 
a,,{x) for xeSniV), i e {1,2}, (23) 
for X e SniT). 

The defined diffusion tensors a^"^ determine differential operators Aq^\x) which approx- 
imate the differential operator Ao{x). Each A'^\x) determines a generalized diffusion 
X (")(•). The generalized diffusions X (")(•) converge in distribution to the original diffusion 
X{-). Therefore it suffices to consider the discretizations yll"-* which are defined in terms 
od a^"'' as described above. The convergence in distribution of diffusions to X{-) 

follows from the convergence of the corresponding semigroups [/'^"•'(•) to the semigroup U{-) 
in C'(K'^). To see this, we first apply Corollarv l2.1l in order to prove the uniform boundedness 
sup(<]^ ||?7("^(t)||*^"-' < /?, where a, /3 do not depend on n. Then we use the standard methods 
of Sobolev spaces in order to prove the convergence of the semigroups in L2{M.'^)- The two 
above properties are combined in order to prove the convergence of semigroups in as 
in the last step of proof of Theorem 14.11 



3.1 Construction for d> 3 

The goal of the overall analysis is to find those discretizations An of the differential operator 
A{x) which have the structure of matrices of positive type. Here we describe a general 
approach, which is based on reduction to a finite number of two-dimensional problems. 

The index set of pairs I{d) = {{ij} : i < j, i,j — 1,2, ... ,d,i ^ j} has the cardinal 
number m{d) = d{d — l)/2. To each index {kl} E I{d) we associate three coefficients, 

Jkl} _ 1 {kl} _ 1 C9^A 

and a two-dimensional bilinear form a^'°'^(-, •), 

a^''^\v,u) = Y I a\f\x)^^v{x)^Ju{x)dx. 

Clearly, for each pair v,u E Cq^''(M'^) the following identity is valid: 

a{v,u) = a^'^'-Hv.u). (25) 

{fe/}e/(d) 

To each of the forms a^^''^{-, •) we associate a sequence of forms ai'^'^(-, •) and matrices aI!'^^ 
constructed by using schemes in two dimensions from the previous subsection. Then the 
matrix 

A^ ^ ^"''^ (26) 

is a discretization of Aq[x). If each aI!'^^ has the structure of a matrix of positive type, 
then An is also a matrix of positive type. However, An can have the structure of a matrix of 
positive type, even though no A^l^^^ is a matrix of positive type. This important property, 
which enables us to construct matrices An of positive type for the case d > 3, can be proved 
from the structure of entries {An)kk±ei- 

First we consider the entry {An)kk+ei defined by (|26l) in the case that < for all 
* 7^ j: *j j = 1; 2, . . . , d. In addition, in order to write down expressions as simply as possible. 
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the index I Cz C is omitted from the notations. The contribution from the sum of entries 
(^An^^^)kk + ei to thc entry {An)kk+ei has the foUowing form: 



— - Oil (^x + ■2(^1 + ^s)^ ~ terms containing a\ 



2, ai3, . . . , aid- 



s>1 



Similarly we can describe the terms containing an for any i — 1, 2, . . . , d, and any / e L. 
The corresponding sum contributing to (A„)fcfc+ei has the following general form: 

1 

ujn{aii,x) -J— J ^ au{hk + hm,i{l,s)), (27) 

where mu{l,s) are defined by the rules of construction of ([M]) . (1^ and (PT|) . The terms 
proportional to 0^5,3 ^ i, are summed with the just defined U!n{aii,x) as shown in the next 
description of the obtained results. 

Discretization procedure 3.1 Let there be given a partition = UiDi into a finite 
number of subsets Di such that all the functions atj are uniformly continuous on each Di, 
and the functions aij,i = 1,2, . . . , d, do not change sign on Di. Let a parameter 

r{l) = {ri(l),r2{l), ■ ■ ■ ,rd{l)) G N'' be assigned to each Di and the matrices An on G„ be 
constructed by the rules i20\) . 1121]) and i26]) . Then their entries have the following properties: 

1. Entries of {An)ki,x — hk, k e Z"^, hi e N{x), are linear combinations ofaij{xij{n, x, I)) 
where Xij{n,x,l) = hk + hmij{l,s), and where the elements mij{l,s) G M!^ for 
i, j, s — 1,2, . . . , d,l Cz L, do not depend on n. 

2. For each grid-knot x — hk: {An)kk ~ — X)j(^n)fcJ- 

3. For each x — hk G cls(£'/) entries in the coordinate directions x ± /le^ are defined by: 



(^")/cfc±e- ~ 73 ^ri{o,ii,x) — ^-jj-\aim{Xim{n,X,l))\ 

7n—l,myti 



4. For each I G L the entries of An in the plane spanned by ei,ej are defined by using 
elements Zij{l) = riiVjBi — rj(l)ej G Z'* (if o-ij < Q on Di) or elements Zijil) = 
r,(0e, + r,(0e, gZ^ (if Oij > on Di), as follows: 

(^")fcfc±..,(0 = - h^r,{l)r,{l) 

An appropriate choice of the parameters r{l) follows from Theorem 13. II 
Some special features regarding the structure of the sets Af{x),x — hk G Gn, should be 
pointed out. If ^ 0,i ^ j, then the maximal number of elements in J\f(x) is 1 + d + d'^. In 
this case the set J\f{x) consists of its center, 2(i-grid-knots in the coordinate directions, and 
2 grid-knots in each of d{d— 1) two-dimensional planes. Since there can be at most two grid- 
knots in a two-dimensional plane (Property 4. of Discretization procedure 13. ip . the entries 
of An'^'^ have the following property. Let the pairs 6^,6^ and es,et span two-dimensional 
planes and let An^'^ , An*'' be the corresponding discretizations which are constructed by 
using parameters r'^''*\ r^**). Then for the construction defined by Discretization procedure 
13. li the following identity (r^'''*))^ — (r^''*^)^ must be valid. 
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3.2 Lower order differential operators 

The discretizations of differential operator B{x) — b{x)W are denoted by The following 
general rule should be obeyed. A positive diagonal entry and a non-positive off diagonal 
entry is associated to each x £ Gn for which b{x) ^ 0. Let us define the sets 3C(i, —) — {x £ 
Gn ■ bi{x) < 0} and analogously X{i,+) — {x & Gn ■ bi{x) > 0}. Then the discretizations 
of {v\Bu) are defined by 

bniv,u) = X! [ X! b,{x)v{x) {a,ih)u){x) + ^ b,{x) v{x) [ai{~h)u){x) . 

These forms have to be summed with the forms (P5|) in order to get discretizations of the 
original form If discretizations (Ao)„ of Ao{x) have the structure of matrices of positive 
type, then obviously {Ao)n + Bn maintain this structure. The so defined discretizations of 
B are usually called upwind schemes. 

The constructed forms a„ of this section are discretizations of the form ([6]). At the 
present level of analysis the constructed discretizations can be justified by the limit a{v, u) = 
lim„ h'^an{v, u), being valid for any pair w, u G Cq^''(M''). 



3.3 Summarized results of the construction 

THEOREM 3.1 Let there be given a partition M"* = U;!?/ into a finite number of connected 
sets Di, each being the union of cubes Gm(jp, x) with some fixed m, so that the functions Uij 
fulfil and the following additional conditions: 

a) The functions Oij are uniformly continuous on Di and aij,i ^ j, do not change sign on 

Di. 

b) For each pair i.j the limit lim|2.|_^oo o,ij{x) has a constant value. 

c) The matrix-valued function x M- d{x), defined by (0), is strictly positive definite on , 

i.e. {z\a{x)z) > j3\z\'^ for some ^ > and all z,x e M.'^ . 

Then there exist discretizations An which are constructed by the rules of Discretization 
procedure \3.U such that An are matrices of positive type. 

Proof: If for each Di we choose the parameters ri{l) so that the entries {An)^i^^^ of 
item 3. of Discretization procedure l3 . 1 l have all negative values, then the rules of construction 
(|2ip ensure the existence of A„ with the structure of matrices of positive type. It remains 
to justify the existence of such parameters r{l) € N"^ for each I S C Let us consider a set 
cls{Di) and the quantity: 



uj(a) ~ inf min min inf inf o,ii(z) 

n leC i=l,2....,d I zeS„{r{l),x}nDi 



nil) ^ (28) 

V ^-TTT sup |a„,(2;)| : x^Gn)- 



zeS„(r(l),x)r\Di 



If uj{a) > 0, the chosen parameters r{l) ensure the positive value of the brackets in item 
3. of Discretization procedure 13.11 If u){a) < 0, the partition should be refined until the 
condition uj{a) > is achieved. In accordance with Lemma l3Tl for each x G G„ there exist 
{ri{x),r2{x), . . . ,rd{x)) such that au{x) - J2m^i{ri{x)/rm{x))\ai,n{x)\ > 0. Due to the 
uniform continuity of the functions Oij on the sets c\s{Di), and b), the described procedure 
results with a desirable result after a finite number of steps. QED 



The uniform continuity of functions Oij on Di, and the inequality w(a) > 0, where io{a) 
is defined by (1281) . imply another important property of matrices An that are constructed 
by the above described procedure. There exists (Tq > 0, independent of n, such that 

> ^, x = hkeGn. (29) 
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The described construction of matrices An for which (j29p holds is called admissible method, 
anticipating that the obtained An have all the necessary properties for the convergence of 
corresponding MJPs to generalized diffusion. Let us recall Ui = Ui{l). 

LEMMA 3.2 Let the matrices An on Gn be discretizations of Aq = —'^dittijdj by an 
admissible method. Then the matrices An are irreducible. If in addition, the matrices An 
are symmetric, then there exist positive numbers M_, M such that the following inequalities 

i i 

are valid uniformly with respect <o n G N. 

Proof: It is easy to check that a symmetric matrix An of positive type is positive 
semidefinite, i.e. (u|^„u) > 0. Let us consider the tensor valued functions a,b where 
a{x) — {aij{x)}fi and b{x) is defined by bij{x) = aij{x) — K6ij. The corresponding auxiliary 
tensors of (U) are denoted by a, b. Due to the strict positive definitness of the matrices a, a 
one can choose k > sufficiently small so that b, b are also positive definite on W^. Let us 
define matrices iJ„ on Gn by the following non-trivial entries: 

i^")kk = 7^' (^'Ofcfcie. = ^7^' « = l,2,...,d. 

We have (u|i/„u) = Yl'i=i \Uiv\2- In accordance with the construction of tensors b,b 
and the inequality (j29p . the matrices Bn — An — nHn are also of positive type. Since the 
symmetric matrix Bn of positive type is necessarily positive semidefinite, i.e. (u|i3„u) > 
for any u G Zo(G„), we have: 

d 

(u|A„u) = K(u|iJ„u) + (u|B„u) > K^|[/iu|^, 

i=l 

proving the lower bound of the assertion. The upper bound follows from (|19p and (|25p . The 
irreducibility follows from the graph theory, since any two grid- knots Xo,y € Gn can be 
connected by a path of the form {xq, Xi, . . . , x^, y} C G„ such that J\f{xk-i) r\J\f{xk) ^ 
for fc = l,2,...,m. QED 

For a differential operator A{x) — Ao{x) + b(a;)V with non-constant coefficients, the 
functions G„ 3 x i-^ S{A,n,x)p{x) of ([T4|) are not identically zero on Gn- In the proof of 
convergence in Section |4] the following weaker result is therefore used: 

LEMMA 3.3 Let the differential operator A{x) satisfy the conditions of Theorem \3. 1\ and 
let the following additional conditions be valid for some a € (0, l); 

a) Functions aij belong to the class (7*^^^"^ on . 

b) Functions bi belong to the class C^"^ on R . 

Then there exists a positive constant k, depending on M, \\b\\ao, Wo-ijW'^^^'^^ and but 
not on n, such that 

sup \S{A,n,x)f{x)\ < k/i" 

xGG„ 

for any f G C(2+")(R'^). 

This result can be easily checked by calculating A„f„ directly. 

Let us mention the following. If the functions aij,bi are uniformly of the class C^^"'""-' 
and C^"^ on Di, respectively, and / £ (7^^+"^ then S{A, n, x)f{x) converges to zero on grid- 
knots X e int(_D/) n ( U„ Gn) ■ Otherwise, n >-)■ 8{A, n, x)f{x) is not bounded as n increases. 
Nevertheless, the convergence in H/j^-spaces |LR2) is ensured as usually. 
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4 CONVERGENCE OF MJPs 



The convergence of MJPs to generalized diffusion is analyzed here in terms of the criterion 
([3]). Therefore, we first need to define explicitly the mappings $„ : Zoo(G„) C{M.'^). 

An element (column) u„ € l{Gn) can be associated to a continuous function on W in 
various ways. In the current setting we define a mapping l{Gn) ^ CiR") in terms of hat 
functions. Let x be the canonical hat function on R, centered at the origin and having the 
support [—1,1]: 

{l + x for X e [-1,0], 
1 — X for X e [0, 1], 
for > 1. 

Then z (j)(h, x, z) = x{h~^iz—x)) is the hat function on M, centered at a; € M with support 
[x — h,x + h]. The functions z i— ipkiz) = OiLi <^(^' Xi,Zi),Xi = hki, i — 1,2, . . . , d, are the 
corresponding d-dimcnsional hat functions with support 5„(l,a;) — J^J + h]. The 

functions tpki') G ^(G„), span a linear space, denoted by E{n,W^). Let u„ e l{Gn) have the 
entries = (u„)fc. Then the function u{n) = X^fce/ ^nkff'k is an element of i?(n,R'') and 
defines an embedding of grid-functions into the space of continuous functions. We denote 
the corresponding mapping by <i>„ : Z(G„) i-^ E{n,U.'^) and write 

u{n) = $„U„ = ^ (u„)^.0;e. (30) 
k 

The inverse mapping : i?(n,M'') i-> ;(G„) is defined by $^^(X]"fc<^fc) = where 
the column u has the entries Uk- It is obvious that the spaces l{Gn) and E{n,R'^) are 
isomorphic with respect to the pair of mappings $„,<I>~-'^. Since hn = 2~", it is clear that 
E{n, R'^) C E{n + 1, R"*) and that the space of functions U„i?(n,R'^) is dense 
[l,oo), as well as in C'(R'^). Let us mention that J^k't'k = 1 on Rf^. For two functions 
v{n),u{n) e E{n,R'^) we have (u(n)|u(n)) = h'^Y.ki-'^kiVkUk where Ski = Wft^kWi^ {4'k\'f>i) ■ 
The numbers Sf-i do not depend on n and the following identity is valid: Ski = 1. Thus 
we have $,7^<i>„ = / in /(G„) and ^n^n^ = / in E{n,R'^). Let P(n) be the projector onto 
E{n,R'^) defined by / P{n)f = $„f„. The mapping can be extended from E{n,R'^) 
to C'(R'^) by defining $"1/ = $-ip(n)/. Thus we have $„$-^ = P(n) in C'(R''). 

If Fn is a matrix on G„, then F(n) — "I>„F„"I>^^ is a linear operator in the linear space 
£'(n,R''). It is easy to verify that ||F(n)||p = for p — l,oo, where the norm || • ||p is 

induced by the restriction of Lp(R'^) to E{n, R''). By applying the interpolation Rietz-Thorin 
theorem [BLj we get ||F(n)||p < |Pri|p for p e [l,oo]. 

The objective of the analysis in this section is the comparison of the Feller semigroup [/(•) 
in G(R'') and the matrix semigroups U{n,t) = $„ exp(-A„i)$-i in E{n,R'^) = E{n,R'^) n 
G(R'^), leading to a proof of (jS)). For this purpose we consider the following initial value 
problems (IVP): 

[dt+ A{x))u{t,x) = 0, u{Q,x) ^ uo{x), , , 

u„(i) + A„u„(t) = 0, u„(0) = Uo„, n e N. ^ ' 

where Uo„ are the discretizations of uq. By using the standard methods in the Sobolev 
space we can first prove a result which is weaker than ([3]). Let A{x) be defined 

by (jlj and Assumption I2.1[ and let An be its discretizations on G„, constructed by the 
rules of Section m We consider u{t) = U{t)f for / e G(R^n W^{R'^), the function f{n) 
defined by (|30p . and u{n,t) = U{n,t)f{n). As proved in [LR2] . for each specified / we have 
lim„_,.oo \\u(t) — u{n,t)\\2,i = 0, uniformly on segments of [0,cx)). 

The function u{t) = U{t)f for / € G(R'^) n WUR'^) is continuous on R"*, uniformly 
on segments of [0,oo), as shown in [LSUj . The corresponding functions u{n,t) are also 
continuous on R , uniformly on segments of [0, oo), as follows from their structure, u{n, t) = 
exp(— A„t)f„. Due to the just described convergence in H^2^(K'^), the sequence il = 
{u{n,t) : n G N} is bounded in W2(R'^), uniformly on segments of [0,oo). Hence, this 
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sequence has a subsequence converging a.e. to u on K'^. However, we need to show the 
uniform convergence. As is usual in such problems, the uniform convergence can be proved 
eventually for an appropriately selected subsequence of il. In our approach, the outline 
of proof of ([3]) is as follows. The original differential operator A{x) is approximated by 
differential operators A(™^(a;) with smooth coefficients. The corresponding semigroups are 
denoted by [/('")(•) in C{M.'^) and C/(™)(n, •) in ioo{Gn)- The limit © is then proved for 
each m G N. Finally, by applying the diagonalization argument to the sequence {^('"^(n, t) : 
m, n € N} we get the desired result. The so outlined steps of the proof are performed in the 
next two subsections. 

Convergence for smooth coefficients 

Let the differential operator Q have the coefficients bi that belong to the class C'-"^ and 
the coefficients that belong to the class C^^"*""-' so that it can be represented as, 

d d 
ij=l i=l 

where b'^{x) = bi{x) — J^j djaij{x). Hence, A{x) can be represented as an elliptic operator 

in non-divergence form with coefficients belonging to the class C'"-* on M'*. This form of 
A{x) makes it possible to use results on the existence of a strongly continuous semigroup in 
(7(2+") (R'') = C(2+")(M^C'(R'*) as developed by Solonnikov [53] (a detailed exposition of 
results can be found in [LSU], Sections 13 and 14 of Chapter 4). Thus we have \\U {t)\\^'^+°''^ < 
exp(crt) for some a >0. 

The matrices An of Section [3] approximate A{x) as described in Lemma 13.31 We have 
($-1^ - A„$-i)/(/ifc) = d{A, n, hk)f{hk) so that: 

\{^-'A~AnK')f\oo<^h^\\f\\^'+"^. (32) 

Now we have a straightforward application of this result on approximations: 

LEMMA 4.1 Lett^U{t) he a strongly continuous semigroup in C'^'^+°'\M.'^) , \\U{t)\\'^+°' < 
exp((Tt) with some ct > 0, such that u{t) ~ U{t)uQ solves the first IVP in i31\) . Let t M> Un{t), 
\Un{t) |oo < 1, &e semigroups generated by —An in Zoo(G„), such that u„(i) ~ J7„(t)uon solve 
the IVPs for ODE in ifg7]j. // the differential operator A{x) on fulfills the condition ifgg)) . 
then the following assertion is valid: For each T > there exists a positive number p{T) 
such that 

sup - U„(i)|oo < piT)\\uof+"h", 

te[o,T] 

for all uo e C^^+"\R'^). 

Proof: The function s H> f(s) = Un{t - s)^n^U{s)uQ e /oo(G„), for < s < 
t,Uo G (7(2+") (M''), has a continuous derivative of the form f'(s) = C/„(t - s) (A„3>~^?7(s) - 
^n^AU{s))un. Therefore the following identity must be valid for each uq G C'-^^"\R'^): 

(^^-^U{t)-Un{t)^-^yQ = j Un{t-s){An^-^-'i>-^A)U{s)uods. 

The /oo(G„)-norm of the integrand is first estimated from above by 

|C/„(t-s)|oo|(A„$-i - $-M)C/(s)uo|oo < k/i"||C/(sK||(2+«) 
< k/i" exp(c7T) ||wo|P+"\ 
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where (|32|) is used. Hence, 



implying the assertion. QED 



Due to the density of C'(^+"^(M'^) in C'(K'') we have the foUowing auxihary result. 

COROLLARY 4.1 Let A{x) be defined by and Assumption \2.1l and let it fulfil the 
conditions of Theorem \3.1\ If the coefficients hi belong to the class C^"-* on R'' and aij 
belong to the class C^^^"'^ on M'' then: 

(i) The operators U{n,t)P(n) converge strongly in C{W^) to U(t), uniformly on segments 

of[0,^). 
(a) The limit is valid. 

The so called classical diffusion, i.e. the process in M.'^ determined by the differential 
operator A{x) = —^ij^iaij{x)didj, is usually simulated by using its representation in 
trems of stochastic differential equations. Corollary 14.11 makes it possible to simulate sample 
paths of classical diffusion in terms of MJPs. This alternative approach to simulation gives 
better results in an estimation of the statistical moments of the first exit time from open 
sets at subsets of the boundary with a rapidly changing normal. 



Convergence in the general case 

Let us define a moUificr x i-> i}{n, x) — h~~'^'d{h~^x) in terms of a non-negative function '&{■) 
of class C^^) on R"^ with the support equal to the unit bah Bi (0) C M'' and / i!){x)dx = 1. A 
smoothing procedure of coefficients a.y , bi is determined by replacing these coefficients with 
the sequence of coefficients a^™' — d{m) * aij,b'f^^ = ^{m) *bi, m <^ N, where * denotes 
the convolution operator. For each m G N the resulting differential operator ^'^'"^(a;) has 
the closure in C'(R'*) and generates the Feller semigroup [/(")(•). Now we consider the 
mappings: 

A{x) I — > {A^"^\x) : m e N}, 
; ; (33) 

An ^ {Ai"^ : m e N}, 

and the operators U^"^^(t) and U^'^\n,t)P{n) in the Banach space C'(]R'*), where the semi- 
group is generated by the closure of ^("'^(a;), and where U^^"\n,t) = exp(— yli^-'t). 
For the double sequence of operators U^"^\n,t)P{n),m,n £ N, we will prove the following 
limits: 

, , 71 — > cx) , , m — > oo 

U'^"'\n,t)P{n) — ^ C/(")(<) U{t), (34) 



uniformly on segments of [0,oo). By using p4p and applying the diagonalization argument 
to the sequence {{/'"'•'(n, ■)P{m) : {m, n) e N^}, we get the main result of this article: 

THEOREM 4.1 Let the differential operator A{x) be defined by Q) and Assumvtion \2. 1\ 
There exists a sequence of pairs [m,n[m)) G N'^ such that 

lim||C/(^^)(i)P(n(m))-C/(<)||oo = 0, 

uniformly on segments of [0,oo). Therefore, the asymptotic ^ is valid for the sequence 
{C^S)(-):meN}. 
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Proof: Since the coefficients 6^™'' belong to the class C*^"-' on M'' and the coefficients 
Oj-™'' belong to the class C^^^""^ on M'', the first limit in (p4| follows from Corollary 14. II 

It remains to prove the second limit in (|34p . that is, the convergence of the semigroups 
[/('"H-) to [/(•) in C{R'^), uniformly on segments of [0,oo). The convergence 
to U{-) in L2(R'^) can easily be obtained (for instance, by using Theorem 6.1, Chapter 1 
in [EKI ) . The convergence of to Uit) in C'(M^), uniformly on segments of [0,oo), 

would follow from such convergence on a dense subspace in C'(M''). We choose the subspace 
&o°'\m.'^) = Co(K'') n C(")(R^) which is dense in both, the space LalR'^) and the space 
C'(M'^). Due to (i) of CorollaryO the operators f/(i) and U'^"'^t) are bounded in the space 
C'(")(M''), uniformly on segments K C [0,cx)), i.e. < ^{K), where P{K) does 

not depend on m. Thus we come to the following conclusion. The operators U(t) — U^"^\t) 
are continuous mappings from C^"^(R'^) into C'(")(M''), \\U{t) - C/(™)(t)||(") < ^{K), and 
they converge to zero in L2(R''), uniformly on segments K C [0, oo). 

Now we apply the following auxiliary resuh to Urn{t) = {U{t) - U^"'\t))v, v G C^°'\r'^). 
Let ii = {un : n e N} be a sequence of continuous functions on [0, 1] x R'' such that: 

a) Unit) G L2(R'^) for each t e [0, 1], and sup{lim„ ||M„(i)||2 : t G [0, 1]} 0. 

b) The functions w„ are uniformly Holder continuous in the following sense. There exist 

a € (0,1) and Cq > 0, which do not depend on t or x, such that the restrictions 
Un{t)\Bi{x) fulfil the following two conditions: 

u4t)\B,{x) e C(")(Bi(a;)), \\un{t)\\^^^ < Co., 

uniformly with respect to a; £ R'^ and t £ [0, 1]. 

Then lim„ M„(t) = in C'(M''), uniformly with respect to t € [0, 1]. 

A proof of this auxiliary result is simple. If the assertion were not valid, there would 
exist a positive number 6 and a sequence of pairs {tk,Xk) G [0,1] x R'^,limfe \xk\ — oo, 
such that Uk{tk,Xk) > S. Due to b) the following must also be valid: \uk{tk)\ > S/2 
on the baU Br{xk) where r = ((5/2ca)-^/". A consequence of these inequalities would be 
\\ukitk)\\2 > 2-^dy/\Br{0)\, contradicting a). QED 

We remind the reader that the approximations A{x) i-^ yl(™)(a;), constructed in the 
proof of Theorem 14. 1[ are not the only ones that are needed for the proof of ([3]) . The 
approximations defined by (1231) are also needed in order to get matrices An of positive type. 



5 SIMULATION OF SAMPLE PATHS 

Here we demonstrate the efficiency of simulation of sample paths of a generalized diffusion 
by using MJPs. We intend to estimate the expectation and the variance of the first exit 
time from a Lipshitz domain. The differential operator A{x) = — 9iaij{x)dj on 

is defined by its diffusion tensor, being a piecewise constant tensor-valued function of the 
form, 

a{x) = plooix), < a^, 

where is a positive number, p is a real number and — (1/4, 3/4)^. Let X(-) be diffusion 
determined by A(x) starting from Xq = (1/2, 1/2). For a bounded Lipshitz domain D the 
expectation and the variance of first exit time from D are given by expressions [Li] : 

£[6*] = Var[6'] = 2||t;||i - £[0]^, (35) 

where u, v are the unique solutions of the following boundary value problems, 

A{x)u{x) — (5(a; — a;o), x e D, A(x)v{x) = u{x), x e I?, ,„„x 
u\dD = 0, v\dD = 0, ^"^^^ 



a{x) 



a a{x) 
a{x) 1 
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and 6{x) is Dirac (5-function at 0. We intend to compute £[6*], Var[6'] by simulations and 
by using deterministic methods formulated in terms of Expressions ((35|) . ([36|l . In order 
to simulate sample paths of X(-) we shall approximate the diffusion by a MJP X„(-) and 
simulate sample paths of In this example we choose — 0.1, p ~ 0.02, D = (0, 1)^ 

and two cases of discretizations, h = 1/200 and h = 1/400. 

The generator of the process Xn{-) in Gn is denoted by A„(gen). Since the entries 
ai2 are non-trivial on Dq, we have to use the construction of Section [3] in order to get 
Anigen) = —An- The parameters ri = 3, r2 = 1 of construction are illustrated in Figure [1] 
These values of parameters ensure A„ to have the structure of a matrix of positive type. 

Two boundary value problems of ([55)) have unique solutions m, u G Li{D) as proved 
in |B01 ILRT| . An efficient numerical method is constructed and the convergence in Li{D) 
is proved in |LR 2j . This numerical method is based on the construction of A„ which is 
described in Section [3l Efficiency of constructed methods is demonstrated by examples in 
which solutions in closed forms are compared with numerical solutions. 

Results of computation are expressed in terms of relative errors: 

<9> 

det 

<9> « » 

det « » sim 

< f >det <SC U »det 

where < >det, < >sim are the estimates of E[0] obtained by using deterministic methods 
(|35|) and Monte Carlo simulations, respectively. Analogously, «C 9 :s>det, <Si ^sim are 
the corresponding quantities for estimates of Var[0]. Some results of computations are given 
in the table bellow. The last column contains the ratios, r = tdet/tsirm of computational 
times tdet and tsim of deterministic and Monte Carlo method, respectively. Sample paths 
are simulated 20000 times. 



h = l/n 


^exp 


^var 


r 


n = 200 
n = 400 


-0.039 
-0.044 


0.018 
0.007 


15.7 
20.3 



Comparison of results obtained by 
deterministic and Monte Carlo methods 

As expected, the first two statistical moments of the first exit time can be estimated by Monte 
Carlo simulations dozen times faster than by using the deterministic method formulated by 
([551) and ([gS]) . 
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